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Abstract 



A new treatment of Bose-Einstein correlations is incorporated in a space-time parton-shower 
model for e + e~ annihilation into hadrons. Two alternative afterburners are discussed, and we 
use a simple calculable model to demonstrate that they reproduce successfully the size of the 
hadron emission region. One of the afterburners is used to calculate two-pion correlations in 
e + e _ —>■ Z° — > hadrons and e + e _ — > W + W~ —>■ hadrons. Results are shown with and without 
resonance decays, for correlations along and transverse to the thrust jet axis in these two classes 
of events. 



1 Klaus Geiger perished in the tragic accident of Swissair flight f f f on Sept. 2, 1998. The three other authors 
dedicate this common work to his memory. 



1 Introduction 



Perturbative parton-shower Monte Carlos [ [T], combined with models for hadronization pro- 
vide a very successful description of experimental data on e + e~ — > Z° — > hadrons, deep-inelastic 
lepton-nucleon scattering, etc.. In most of the applications made so far, attention has been con- 
centrated on distributions and correlations in momentum space. However, there are some key 
aspects of the physics where better understanding [ |3| of the space-time development of the 
hadronic system is desirable [ |]]. This is particularly true for the treatment of dense hadronic 
media, such as those produced in heavy-ion collisions, where the formation and expansion of 
the system are of both experimental and theoretical interest. A prototype for the treatment of 
such questions may be provided by the reaction e + e~ — > W + W~ — ► hadrons, where the W ± 
do not decay independently, but in a hadronic environment created by each other. This may 
engender collective effects such as colour reconnection [ |5], [7], ^| , parton exogamy [ || [TTJ and 
Bose-Einstein correlations [ [11], |12|> 0] that ma y be detected by experiment [ [L4], [15], |16| , and 



could be of relevance to the measurement of mw at LEP 2 [ |17|] . 

A parton-shower Monte Carlo has recently been developed [ p| which incorporates the 



information on the space-time development that is encoded in perturbative QCD [ 19], and 
combines it with a phenomenological spatial criterion for confinement [ [20j to provide a com- 
plete space-time description of hadronization. This tool has been applied to the analysis of 
e + e~ — > Z° — > hadrons [ [|, e + e~ — > W + W~ — > hadrons [ || fLO] , deep-inelastic lepton- 
nucleon scattering [ |2T| and relativistic heavy- ion collisions [ |22|, 1231, UM. In the application 



to e + e~ — > W + W" — > hadrons, it has provided new insight into collective effects such as 
parton 'exogamy' [ [TTJ, namely the marriage of partons from different parents to produce 
daughter clusters of final-state hadrons. In the application to relativistic heavy-ion collisions, 
it has provided useful insights into such issues as the formation and local thermalization of 
the dense nuclear fireball, hadron production [ |2Jj] , and the possible suppression of the J / ip [ 



23 1 . However, little attempt has so far been made to incorporate Bose-Einstein correlations 



into this space-time model in a realistic way. 

Bose-Einstein correlations have been analyzed in many experimental situations, including 
e + e~ annihilation [ [14], [15|] where there has also been considerable recent theoretical progress [ 



2JJ, and have been used extensively as a tool to analyze the hadronic fireballs produced in 
relativistic heavy- ion collisions [ |26], |27j. Considerable recent progress has been made in the 
development of the formalism for analyzing Bose-Einstein correlations [ and for implement- 



ing them in an afterburner for models of hadron production [ 120, BUI [31]]. It was shown that 



Bose-Einstein correlations in the two-particle momentum spectra allow for a detailed recon- 



struction [ ^7], of the geometry and dynamical state of the reaction zone from which the 
final-state hadrons are emitted. The purpose of this paper is to describe the implementation 
of Bose-Einstein correlations in the space-time parton-shower Monte Carlo mentioned above 
(see also [ 0), and to describe pilot applications to the reactions e + e~ — >• Z° — > hadrons and 
e + e~ — ► W + W~ — > hadrons. This work should pave the way for a detailed space-time analysis 
of hadron production in these reactions using data on two-particle momentum correlations. 

We introduce in Section 2 of this paper "classical" and "quantum" algorithms suitable for 
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the calculation of Bose-Einstein correlations in Monte Carlo codes for hadron production. Both 
algorithms differ in how the numerical event simulation is used to define a quantum mechanical 
phase space density of emission points. We test these algorithms in Section 3, using a simple and 
analytically solvable model [ Q for the hadron source. We verify that Bose-Einstein analysis 
tools applied to the hadronic spectra generated by the two versions of the afterburner reproduce 
correctly the input source geometry. Then, in Section 4 we apply the "quantum" afterburner to 
hadronic Z° and W + W~ final states generated by a parton-shower Monte Carlo. We calculate 
two-pion correlations with and without resonance decays, studying both the longitudinal and 
transverse momentum dependences of the correlation functions. 

In a final Section, we mention possible future studies using the approach introduced in this 
paper. These would include implementation of the "classical" version of the Bose-Einstein 
afterburner, and exploration of the influence on the "quantum" afterburner results of varying 
the assumed wave-packet size. 



2 Bose-Einstein algorithms 

In this section we explain the algorithms with which we later calculate two-particle correlations 
of identical pions from perturbative parton-shower Monte Carlos. 



2.1 General considerations 

Bose-Einstein correlations reflect the phase-space density of the hadronic source created in the 
collision. Contrary to single-particle momentum spectra, they thus also provide access to the 
space-time structure of the reaction zone. A consistent numerical simulation of Bose-Einstein 
effects on the two-particle and many-particle momentum distributions thus requires by necessity 
an algorithm which propagates the particles in phase space, rather than in momentum space 
only. This is what the parton-shower cascade event generator VNI does. 

The Bose-Einstein symmetrization effects result in an enhancement at small relative mo- 
menta q of the 2-particle coincidence spectrum relative to the product of single particle spec- 
tra. The q range of this enhancement is inversely related to the size of the emission region 
in space-time. All existing shower Monte Carlos, whether formulated in phase space or only 
in momentum space, are based on a probabilistic description and thus do not correctly de- 
scribe the many-particle symmetrization effects of the quantum-mechanical time evolution. 
The corresponding quantum-statistical corrections must therefore be implemented, in some 
approximation, by an "afterburner" at the end of the classical time evolution. 

In momentum-space based Monte Carlos like JETSET, one tries to implement the clustering 
of identical bosons at small relative momenta by shifting the final state momenta according to 



certain prescriptions [ O, O. These shifting prescriptions are not unique and lead to changes 



in the invariant mass of the particle pair, and thus do not allow one to conserve simultaneously 
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energy and momentum. Furthermore, they involve a weighting function which is put in by hand 
but should, in principle, reflect the (unknown) space-time structure of the simulated event. A 
recent attempt to relate the weighting function directly to previously unused information on 
the space-time structure of the particle-production process is described in [ However, 
its connection to the position of the hadrons at "freeze-out", i.e., decoupling from the strong 
interactions, remains at most indirect. 

In the present paper, we study two algorithms [ [ftj to implement Bose-Einstein corre- 
lations at the end of the Monte Carlo simulation. These algorithms do not shift the particle 
momenta, nor do they alter the output of the event generator in any other way. They calcu- 
late the single-particle inclusive momentum distribution directly from the output momenta of 
the generator, and the two-particle coincidence spectra from the space-time coordinates and 
momenta of particle pairs from the generator output. They differ in the way in which they 
associate with the event generator output a quantum- mechanical Wigner phase space density 
S(x, K). Both algorithms assume that the particles propagate freely from the source to the de- 
tector and include only the quantum-statistical pairwise correlations between identical bosons. 
Generalizations of these algorithms to include final-state interactions [ |35| and multiparticle 
correlation effects [ ^] have been proposed but not yet implemented numerically. 

The two-particle correlation function is constructed as the ratio of the two-particle coinci- 
dence spectrum P2(p a ,p b ) and the product of single-particle inclusive spectra, P\(p ab ), 

where q = p a — p b is the relative and K = (p a + p b ) /2 is the average pair momentum. With 
the assumption of independent particle emission the two-particle correlation function (HD can 
then be written as [ [J7|, [38|, |3£ 



C(q, K) - Af s \1 + s ^ ^ ^ s ^ j , (2) 

where S(x,p) is the single-particle Wigner phase-space density of the source. In this work we 
choose the normalization M = J\T S = 1 in presenting our results. The implications of other 
choices of normalization are discussed in Section |2.2| . The four-vectors p a ^ b in the denominator 
on the r.h.s. are on-shell while the numerator contains the off-shell four-vectors q and K with 
q° = E a — E b and K° = (E a + E b )/2. The main question is how to relate the event generator 
output to this Wigner density, and how to simulate the r.h.s. of Eq. (0). This will be discussed 



in Section 2.3. 



2.2 Normalization of the correlator 

The normalization M in ([I]) does not affect the space-time interpretation of the correlator, and 
the reader who is only interested in the latter can skip the present subsection. The subtle point 
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we discuss here is that, in the context of event generator studies, the normalization TV of the 
correlator is only fixed after requiring that the Bose-Einstein algorithm affects the simulated 
multiplicity in a particular way We start by recalling the quantum field-theoretical definitions 
of the single- and two-particle spectra, 

Pi(p) = E p (a%), (3) 
Pi(Pa>Pb) = E aE b (d p a p a Pb a Pa ) , (4) 

where (...) indicates the ensemble of physical states (events) for which the correlator is calcu- 
lated. This implies the normalizations 

J p iiP)^ = (N), (5) 
fP2( Pa ,P b )^^ = (iV(iV-l)), (6) 

where N = f(d 3 p/E)a p a p is the particle number operator. We now discuss the physical 
implications of two different normalizations of the correlator: 



1. One can interpret the two-particle correlator as a factor [ [13 



d 6 a™/d 3 Pl d 3 p 2 = C(q, K) d«a™/d 3 Pl d 3 p 2 (7) 

relating the measured two-particle differential cross sections on the l.h.s. to the differential 
two-particle cross section resulting from the simulation. Requiring that the Bose-Einstein 
afterburner conserves event multiplicities on an event-by-event level, the corresponding 
momentum integrated total two-particle cross sections have to coincide, a^f = u^®- This 
is the appropriate starting point if total pair cross sections are used in the tuning of the 
event generator which then, of course, should not be changed by the Bose-Einstein after- 
burner. The normalization satisfying these requirements normalizes both the numerator 
and denominator of ([[]) separately to unity [ , 

M - 



(JV(iV-l)} 

This results in a normalization J\f s < 1 of the two-particle correlator (J|) [ [41], |3(| 42 |. 
2. A different choice of normalization often used in heavy-ion physics is [ 

M = M s = l . (9) 

Combining ([TJ) and (0), it follows that 

d 6 N I d 3 N\ ( d 3 N\ 

ftfe., ft) = E. E il ^ i > P l( pJ P l(Pi ) ={E« W \ [E iWi , (10) 
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and, because of (|5]||), also that 



{N{N- 1)) > (iV) s 



11) 



If we interpret the r.h.s. of these equations as the pair spectra and pair multiplicity from 
the event generator, implying that the generated multiplicity has a Poisson distribution, 
((N(N — l)) gcn = (N)g en , then this implies that the Bose-Einstein effects have increased 
the pair multiplicity. This may account for some of the effects of Bose-Einstein statistics 
on the particle-production processes prior to freeze-out [ H3 



Depending whether we require for the Bose-Einstein afterburner the conservation of event 
multiplicities on an event-by-event level, or aim to mimic Bose-Einstein effects during the 
particle-production processes as well, the normalization of the two-particle correlator is thus 
either smaller than unity or unity itself. In the present paper, we are only investigating the 
space-time interpretation of the two-particle correlator, and hence we can set M = J\T S = 1 
without any loss of generality. 



2.3 Wigner densities and event generator output 

We now explain how we construct a two-particle spectrum with the properties fl9[-|TT|) from the 
event generator output. For simplicity, we discuss only one particle species, say n + . The event 
generator yields for each collision event m a set of final (on-shell) 7r + momenta pi = (Ei,Pj) 
and last interaction points = (ti, r^), with i = 1, 2, . . . , N m where N m is the total number of 
7r + created in event m: 

{{r i ,p i )\i= 1,2,...,AU- (12) 
They define a classical (positive definite) phase-space density 

1 iVevt 1 iVcvt N m 

PoUx,p) = — £ P £L(x,p) = r EE ^ (4) (* - n) 5 {4 \p - Pi ) • (13) 

JV evt m= i iv evt m= i j = i 

The distributions p^l s (x,p) for individual events cannot be taken as Wigner densities since 
they fix the particle coordinates and momenta simultaneously, thereby violating the uncer- 
tainty relation. This can affect the calculation of the two-particle correlator significantly [ |31 



Furthermore, p c \ass(x,p) is always positive, whilst the Wigner density S(x,p) can, at least in 
principle, become negative. Only when averaged over sufficiently large phase-space regions 
is the latter guaranteed to be positive definite. On the other hand, it is unlikely that such 
Zitterbewegung oscillations of S(x,p) or the spiky structure of p c iass(. x ,p) affect the correlation 
function at small q where the Bose-Einstein effects become visible. It is well known [ |2S] that 
the width of the correlation function reflects only the r.m.s. width of the Wigner density S(x,p) 
in coordinate space, and that finer structures in S(x,p) (like spikes or quantum oscillations) 
show up in the correlator only at large q and are very hard to resolve experimentally. Further- 
more, since the event generator performs a Monte Carlo simulation of a dynamical evolution 
which is based on quantum-mechanical transition amplitudes, averaging its output over many 
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simulated events should generate a smooth phase-space distribution (|T^) which is not in conflict 
with the uncertainty relation. 

Following these arguments, one can try to identify directly the classical phase-space density 
Pciass(^,p) fll3D , averaged over sufficiently many events, with the on-shell source Wigner density 
S(x,p) in (Q), in the following sense: 

Pciass(x,p) = 26(p°) 5(p 2 - m 2 )S(x,p) . (14) 

This ensures the correct normalization to the average multiplicity (N): 

[fl fd 4 xS(x;p,E p ) = [d 4 pd 4 xp dass (x,P) = T^i2N m =(N). (15) 

J rjp J J iVevt m=1 

The identification fllil) gives rise to the "classical" version of our Bose-Einstein afterburner [ 



2J|, to be discussed in Sec. |2.3.1 



Alternatively, if one wants to avoid the conceptual difficulty of relating an expression like 
fljjD, where every term under the sum explicitly violates the uncertainty relation, with the 
source Wigner density, one can associate the set of phase space points ( |T2D with the phase- 
space locations of the centers of minimum-uncertainty wave packets [ |30| : 

(r i>ftl tO — fox,*) = —^—e-^ 2 ^^ x . (16) 

In this case one enforces quantum-mechanical consistency of the emission function S(x,p) at the 
level of each individual simulated event. The identification ([if]) gives rise to the "quantum" 



version of our Bose-Einstein afterburner [ ^9], [JIJ, to be discussed in Sec. |2.3.2| . The word 



"quantum" in this case stresses the quantum-mechanical consistency of the Wigner density 
on the event-by- event level (which may indeed be requiring too much), while the "classical" 
algorithm generates a quantum-mechanically consistent emission function only on the ensemble 
level, and only if p c i ass does not violate the uncertainty relation (see Sec. §). 

Before turning to a discussion of these two algorithms we shortly comment on the underlying 
assumptions. The use of single-particle Wigner densities S(x,p) implies that the iV-particle 
production amplitude factorizes into one-particle production amplitudes [ f|4], |28| . In general, 



P^iPaiPb) ls given by a sum over the two possible permutations of the Fourier transform of the 
quantum mechanical two-particle Wigner density S2(x a ,p a ; Xb,Pb) of the source at freeze-out [ 
PI ; here we assume S2(x a ,p a ; Xb,Pb) = S(x a ,p a ) S(xb,Pb)- This assumption (which amounts to 
a Wick decomposition of the r.h.s. of [ p8|| ) implies that the two particles in the pair are 
emitted independently from each other. It thus neglects dynamical correlations between the two 
particles in the pair, due, e.g., to energy-momentum conservation, as well as certain quantum- 
statistical correlations which may be induced on the two-particle level by the symmetry of the 
multi-particle final-state wave function. While the neglect of dynamical correlations is probably 
well justified for heavy-ion collisions for which our algorithms were developed [ ^9], |30|, the same 



is much less obvious for e + e collisions. At high energies, however, we expect such dynamical 
correlations to affect the two-particle spectrum mostly at large values of q, where kinematical 
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constraints play an important role, and not to interfere with the Bose-Einstein correlations at 
small q. If this is true, they cancel from the ratio ([[[) as constructed by our algorithm (see 
below). Multi-particle symmetrization effects, on the other hand, are more of an issue in heavy- 
ion physics [|3(| f|6| where the rapidity densities of the produced particles are large, while their 
neglect in e + e~ collisions seems unproblematic. Furthermore, it is known [ |42| that for certain 
classes of multiplicity distributions they do not destroy the factorization of the two-particle 
Wigner density which is assumed here. 



2.3.1 "Classical" version of the Bose-Einstein afterburner 



We start from (|13l) and fll4|) . The momenta returned from the event generator are on-shell, and 
we hence write from now on S(x, p) respectively p c \ass{x,p) for the on-shell distributions. The 
5-function structure of p c i aS s requires one in practice to bin in the momentum variable (since 
the x-dependence is integrated over in (0), no binning in x is necessary there). For this purpose 
we introduce the normalized "bin functions" with bin width e 

S (e) = f 1/e 3 : Pj ~ | < Pi,j <Pj + ^ {j = x, y, z) 
Pi ' p 1 : otherwise 

or, alternatively, properly normalized Gaussians of width e, 

^ = -4p ex p(-(ft-p) 2 A 2 l • fi 



In the limit e — > 0, these Gaussian bin functions reduce to the properly normalized 5 functions 
$ (Pi ~ P)- For each event m we calculate the numerator and denominator of (|7|) separately. 
We find for the invariant two-particle spectrum in the numerator [ ^] 

N m 

P2( Pa , Pb ) = E a E b J2 {fi£, Pa fi£l Pb + fi^K fiiU cos(g • (r 4 - r s ))) (19) 
and for the product of single-particle spectra in the denominator 

Pi (p a ) Pi (p b ) =wEE fi£, Pa fil% ■ (20) 

In (|T~9|) , K=(p a +p b )/2 and q=p a —p b define the point in momentum space at which the corre- 
lator is to be evaluated. Please note that the momenta pij of the generated particles determine 
only which pairs are selected and contribute to the correlator, but their weight in the correlator 
(in particular the cosine in the exchange term) depends only on the space-time coordinates and 
not on the momenta of the generated particles. 

The correlator ([[]) is obtained by averaging the numerator and denominator separately 
over all events, J2m=i ■ ■ ■■> anc ^ t nen taking the ratio. Direct insertion of ( |T3| , |n| ) into 



gives fll9l , |20D without the restriction i ^ j on the summation indices. This is a discretization 
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artefact, and the pairs with i = j formed from the same particle must be removed by hand in 
this approach. To preserve the normalization of the correlator we also remove them from the 
denominator. Replacing cos(g ■ (r, — r,)) by exp(iq ■ (r$ — rj)), which is allowed by symmetry 
under the exchange % «-> j, the weight function can then be factored, and we obtain 



C(q,K) 









2 




2" 






7WV m r(e) > 

ylui=l °P l ,PaJ 




\^N m r(e) 

2^i=i °Pi, 


Pa U Pz,Pb 



(21) 



The subtracted terms in the numerator and denominator remove the spurious contributions 
from pairs constructed of the same particles. The factorization of the weight function provides 
a dramatic simplification. Each of the sums in (|2T| ) requires only 0(N m ) manipulations, a 
clear advantage for large average event multiplicities (N) over the evaluation of (|19D , which 
involves O(N^) numerical manipulations. Unfortunately this fails once final-state interactions 
are included, since the corresponding generalized weights no longer factorize [ |35| . Also, if 
one wants to account for multiparticle symmetrization effects, more than 0{N^ n ) numerical 



manipulations are typically required [ |36 



In general the result for the correlator at a fixed point (q, K) will depend on the bin width 
e. Finite event statistics puts a lower practical limit on e. In practice the convergence of the 
results must be tested numerically. We discuss these statistical requirements in Sec. |3] in detail 
for a toy model. 



2.3.2 "Quantum" version of the Bose-Einstein afterburner 



In the "quantum" version of the Bose-Einstein afterburner, the phase-space coordinates (ij, r iy p t ] 
of the generator output are interpreted as the centers of normalized minimum-uncertainty Gaus- 
sian wavepackets (IIB) of spatial width a. For the one- and two-particle correlator, one finds 
instead of (EB [ m m 



E, 



d 3 N 

C{q,K) 
Siip) 



rp Ncvt rp N cv t Nm 



N. 



evt m= i 



evt m= i j = i 



1 + e 



-*V/2 



E 



iVevt 

m=l 



Em=l "m(Pa) "m(p b ) ~ Ei=l Si(p b ) 



(22) 

(23) 
(24) 



This result can be derived either directly from ([16]) following [ |3(| , or by replacing the products 
of S functions in (|i~3"D by the Wigner densities of the corresponding wave packets, identifying the 
Wigner density S(x,p) as the sum of the corresponding individual Wigner densities [ ^9], pOf] : 



pS1(*,p) ^ S^ m \x,p) = ^£5(*° -U) exp 



1=1 



(25) 
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In the second derivation, based on (]25f) , the spurious contributions from identical pairs must 
again be removed by hand. Now (p5|) is correctly normalized to the number of particles N m in 
the event: 

jfl Jd i xS^\x,p) = N m . (26) 

The Gaussian single particle probability Si(p) describes the contribution of the generated par- 
ticle i to the momentum spectrum at p. In the quantum algorithm, it is the counterpart of the 
bin function dL e ) p in the "classical" afterburner. The limit of vanishing bin width e — ► cor- 
responds to the limit a — > oo in which the wavefunctions (|16"D become momentum eigenstates. 
The difference between the two algorithms is then essentially the prefactor exp(— a 2 q 2 /2) in 



fl23|) which is a genuine quantum contribution. A momentum eigenstate is infinitely delocalized 
in space, and the prefactor exp(— a 2 q 2 /2) ensures that this infinite source size is reflected in 
a sharp correlator C(q, K) = 1 + 5 q ,o- We emphasize that while e — > is the relevant phys- 
ical limit for the "classical" afterburner, a — > oo is not the relevant limit for the "quantum" 



algorithm (|22|)-(|2 



It might seem natural to interpret a in terms of the size of the hadronic cluster at its 
formation, or of its wave function at decoupling from the other particles, which would suggest 
a values in the range a ~ 1 fm. However, such arguments are not rigorous, and in the present 
paper we treat a as a phenomenological parameter. One could in principle adjust this parameter 
as part of a general optimization or tuning of the Monte Carlo, but such a study extends beyond 
the scope of this paper. 

It is important to note that, due to the smooth intrinsic momentum dependence of the 
Gaussian wavepackets (|16D, the correlator ( f2~3|) is a continuous function of both q and K, even 
though the event generator output is discrete. On the other hand, due to the piecewise constant 
nature of the bin functions (|17])) the correlator ([H]) is only a piecewise constant function of its 
arguments, which may, in practice, require binning in both q and K. 



The Bose-Einstein algorithms explained in section |2.3| can be applied to any model that 
gives a particle phase-space distribution, irrespective of the dynamical history of the particles. 
The aim is to reconstruct from the Bose-Einstein correlations information about the space- 
time history of the dynamical evolution, as one attempts to reconstruct in real-life experiments 
the space-time structure of collisions from the particle distributions measured by a detector. 
However, with a particle sample from an event generator model, detailed knowledge about the 
dynamical evolution is available. This allows one to cross-check whether the generated dynamics 
reproduces the measured Bose-Einstein correlations, i.e., this provides an experimental test of 
the generated space-time interpretation [ B3]. 



3 Tests of the Bose-Einstein afterburners 



In this section we show numerical tests of our afterburners using a simple toy model for the 
source which allows for analytical calculations of the correlation function. We thus illustrate 
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the algorithms discussed in section 2J3 before turning in Sec. [| to realistic parton-shower cal- 
culations. 



3.1 Analytical model studies 



We explore the above algorithms with a simple model emission function first proposed by Zajc 



PcZ s ( x >P) = A/; exp 



1 / x 2 x • p p 2 



- 2s— ^ + 



2(1 -S 2 )\R 2 i2o p Q p2 



S(x°) , < s < 1, (27) 



N 

K = E v (28) 

This distribution is normalized to an event multiplicity N, and is localized within a total 
phase-space volume 

V p . s . = (2R S P ) 3 R s = RoVT^ 2 , (29) 

which vanishes for s — > 1. This s-dependence allows one to study the performance of our 
numerical algorithms for different phase-space volumes. The parameter s smoothly interpolates 
between completely uncorrelated and completely position-momentum correlated sources: for 
s — > 0, the position- momentum correlation in ( p7[) vanishes, and we are left with two decoupled 
Gaussians in position and momentum space. In the opposite limit the position-momentum 
correlation is perfect, 

limpS^P) ~ 5 {3) - 5(x°) , (30) 

and the phase-space localization described by the model violates the Heisenberg uncertainty 
relation. 

How are these properties reflected in the one-particle spectra and two-particle correlation 
functions? It turns out that in the Zajc model the two-particle correlator is independent of 
the pair momentum K, irrespective of s. Due to the spherical symmetry of the source and its 
instantaneous time structure, the correlator is thus characterized by a single, l^-independent 
Hanbury- Brown- Twiss (HBT) radius parameter. 

In the "classical" interpretation S(x,p) = p^^ s (x,p), and the one-particle spectrum and 
two-particle correlator read 

r, dN n N ( P 2 \ 

E ^ P = Ep j2^W^ eXP \^J ' (31) 
C(q,K) = l + exp(-R 2 class q 2 ) , (32) 

^lass = Rl (l - j2R~Pof) ■ (33) 
For sufficiently large s, when the phase-space volume becomes smaller than unity, 



,> '- = f-W** ll ' <1 ' (34) 
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the HBT radius parameter turns negative, which leads to an unphysical rise of the correlation 
function with increasing q 2 . Since K,. s . = 1 corresponds to the volume of an elementary phase- 



space cell, the change of sign in ( p3|) is directly related to the violation of the uncertainty 
relation by the emission function (|27|). 



In the "quantum" interpretation, phials gi yes the distribution of centers of Gaussian wavepack- 
ets, and the Wigner phase space density is obtained from phials Yia ° (p6|). The one-particle 
spectrum and two-particle correlator then read 

dN N ( p 2 \ 

C(q,K) = l + exp{-q 2 R 2 qm \ , (36) 

1 



Rim - R 2 1 - , 9 ppN 2 > ( 37 ) 



R 2 = R 2 + y> p2 = p o+^- ^ 

In this case, R and P satisfy 2RP > 1 independent of the value of a, and the radius parameter 
R 2 m is now always positive. Even if the classical distribution Pciass^f) violates the uncertainty 
relation, its folding with minimum-uncertainty wave packets leads to a quantum-mechanically 
allowed emission function S(x,p), and to a correlator with a realistic fall-off with q 2 . The 
limiting cases are also as expected: For R — > oo the source extends to spatial infinity and the 
correlator collapses to a Kronecker 5 function at |q|=0. For Pq — > oo, the source is momentum- 
independent, and the HBT radius measures a combination of the geometric extension of p c iass 
and the spatial wave-packet width a: R 2 ^ = R 2 + cr 2 /2. For s = 0, one recovers the expressions 
given in [ [30]]. The folding with wavepackets modifies the geometric size of the source by adding 



in quadrature the intrinsic width of the wavepacket, R 2 ntr =a 2 /2 and the size of the classical 
distribution p c i ass . The extra term is exactly reflected by the prefactor e _CT q I 2 by which ( |22| ) 
differs in structure from the classical result (|22| 



However, the spread of the one-particle momentum spectrum (^) receives an additional 
contribution l/2cr 2 . Choosing a too small increases this term beyond phenomenologically rea- 
sonable values, whilst choosing it too large widens the corresponding HBT radius parameters 
significantly. This restricts the range of phenomenologically acceptable a values for the "quan- 
tum" version of the Bose-Einstein afterburner. 



3.2 Event generator studies 

To test the afterburner, we have mimicked the role of an event generator by creating a Monte- 
Carlo phase-space distribution of iV phase-space points {(ri,p i ,ti)\i = 1, . . . ,N} according to 
the distribution p^als i n (FD - This Gaussian model distribution allows one to compare the 
numerical results of the Bose-Einstein algorithms to the analytical expressions obtained above, 
thus testing statistical requirements, the accuracy of the numerical prescriptions, and the role 
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of the bin width in the "classical" algorithm. Its generic properties in both the "classical" and 
"quantum" versions can be read off from Fig. [1|. 




Figure 1: Generic properties of the one- dimensional Zajc model, a): The HBT radius (|33|) from 
the "classical" algorithm as a function of the strength s of the position-momentum correlations 
in the source. Different curves correspond to different combinations of the model parameters 
Ro and Pq and to different bin sizes e. b) : s-dependence of the HBT radius parameter ( |37| ) for 
the "quantum" algorithm. Different curves are for different combinations of model parameters 
Rq and Pq and different wave packet widths a. c) and d): The two-particle correlator in the 
"classical" and "quantum" versions of the afterburner, for different sets of model parameters. 
The numerical results were obtained by analyzing 50 events of multiplicity 1000, i.e. iV" pa i rs = 
2.5 * 10 7 . They show small deviations for the "classical" version, but coincide within the line 
widths for the "quantum" version. 

The HBT radius parameter (|33p of the "classical" prescription, depicted in Fig. |I]a, strongly 
depends on the position- momentum correlation in the source. Below s crit , it is positive, which 
corresponds to a quantum-mechanically allowed Wigner function. Above s crit , -Rciass turns 
negative, i.e., -R c iass is imaginary. The value of s cr i t depends on the total source size 2RqP in 
phase space. In Fig. |l]a we exploited this by varying P between 80 and 200 MeV, keeping 
Ro = 1 fm fixed. In the plot one sees again that the HBT radius parameter takes unphysical 
imaginary values as soon as the phase-space volume (2R S P Q ) S becomes smaller than 1. 



As explained in Sec. [2.3. 1| above, the "classical" Bose-Einstein algorithm requires a smearing 
of the momentum-space S functions in ( jT3]) by bin functions ( [T7| ) or ( jT8]) of width e. The physical 
situation is recovered in the limit e — > 0, but a careful investigation of this limit is numerically 
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difficult. However, for the Gaussian bin functions (18) the HBT radius parameter can be 
obtained analytically for finite bin width e: 

^ciass(e) = I ~ e 2/(2p2) ( X + 2P 2 (1 - S 2) ~ (2i? s P ) 2 ) • (39) 

Comparison with fl33"D shows that the numerical results should be close to the physical ones if 
one chooses 

e<^V2P . (40) 

This provides the useful information that in practice the scale for e is set by the width Pq of 
the generated momentum distribution, independent of the geometric source size R s . 

In Fig. |l|a we have also plotted the e dependence of the HBT radius parameter. Clearly, 
for fixed bin width the approximation of the true HBT radius parameter (|33|) becomes better 



with increasing Pq, as suggested by (39). More generally, the net effect of a finite bin width is 



always to increase the apparent size of the source. 

In Fig. |l]b we show the HBT radius obtained from the "quantum" version of the afterburner. 
Now the situation is qualitatively different: the HBT radius is always positive, since the smear- 
ing with Gaussian wave packets always ensures consistency with the uncertainty relation, and 
its s-dependence is much weaker since the wave packets smear out the unphysically strong 
position-momentum correlations in p^^ s - The different curves shown in Fig. |l]b illustrate, for 
fixed classical source radius Rq, the dependence of the HBT radius parameter on the width P 
of the classical momentum distribution and on the wave packet width a. Wave packet widths 
o > Ro not only change the HBT radius itself, but also its dependence on P significantly. 

In Figs. |l]c,d we present for characteristic model parameters the corresponding two-particle 
correlation functions. The analytical curves (|32| ) and fl36|) are compared to numerical results 
from the algorithms (|2lD and (^) applied to a Monte-Carlo distribution of phase-space points 
{(?*», p 4 , ti)\z = 1,...,N} obtained from the distribution phials * n (ffiD- The plot shown used 
N evt = 50 events of multiplicity N m = 1000. We emphasize that only the total number of pairs 
in the event sample, ^N evt N m (N m — 1), is statistically relevant. Our choice of N evt and N m 
hence illustrates the properties of the algorithms for both high and low multiplicity events. 

For the "quantum" algorithm, the numerically simulated correlator in Fig. [Tjcl coincides 
with the analytically calculated one (|36|) within the line width. Small differences between 
the analytic and numerical results are seen for the "classical" algorithm in Fig. |l]c. In order 
to understand these differences in the performance of the two algorithms quantitatively, we 
have studied their statistical requirements in the following way: from the distribution phials 
we generated a set of iV fit = 5000 samples of N evt events, each event containing N m = 100 
particles. For each of the 5000 event samples, we calculated the two-particle correlator with 
both algorithms and determined the HBT radius R^ from a Gaussian fit to this correlator. The 
statistical deviations from the analytically known exact results -R c i aS s(e) and R qm , respectively, 
were then determined as a function of the sample size iV evt via 

1 Nat 

A stat (iV evt ) = — E ( R M - P exact ) 2 • (41) 
n=l 



iVfit 
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In general, increasing the event multiplicity N m or the number of events N evt improves the 
performance of the algorithms. Here we focus on the typical situation that the average event 
multiplicity N m is fixed by the simulated physics, while the number of events in the event sample 
can be increased by a longer running time of the (numerical) experiment. The corresponding 
statistical performance of both algorithms, measured in terms of A sta t(-Wevt), is plotted in Fig. ||. 
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Figure 2: The average statistical deviations A stat (iV evt ) of the "classical" and "quantum" 
Bose-Einstein algorithms as a function of the number of events in the sample. Here, the event 
multiplicity is N m = 100, and only the total number of pairs per event sample is statistically 
relevant. 

For the "quantum" algorithm, the statistical fluctuations decrease like A stat (iV evt ) ~ -J— . 

iVevt 

Also, their absolute value is small: for only N evt = 10 events, the fluctuations in the fitted 
values .Rfit are already smaller than 0.1 %. This is the reason why in Fig. |T|d for the "quantum" 
algorithm the simulated values coincide so well with the analytical ones. We also observe that 
A s tat(^Vevt) increases for larger values of a. The reason is that A sta t(^evt) oc R 2 x&ct) which 
increases significantly with increasing a (see Fig [l|b). The normalized fluctuation measure 
A sta t(^evt)/-Rexact decreases slightly with increasing cr, since the finite wave-packet width smears 
out the discrete classical emission function flT5| ) and thereby reduces the statistical fluctuations 
in the algorithm. 

In comparison, the "classical" algorithm shows statistical fluctuations which are approxi- 
mately two orders of magnitude larger. One sees clearly how A gtat (iV evt ) increases, i.e., the 
statistical requirements increase, if one goes to smaller bin widths e, as needed to realize the 
physical limit e — > 0. Also, at least for small values of iV ev t < 100, the fluctuations A s t a t(A r ev t) 
decrease more slowly than 1/N evt . There are several reasons for these differences between the 
"classical" and "quantum" algorithms. Numerically, we observe that, in the "classical" algo- 
rithm, the simulated correlator ( [21] ) has even for the present Gaussian model a tendency to 
become non-Gaussian. This is seen, e.g., in the slight deviations in Fig. |l|c for s = 0.2. These 
non-Gaussian effects depend on N evt and manifest themselves in the slight wiggle in Fig. Q in 
the curve corresponding to e = 40 MeV, which is a relatively large bin width. Secondly, we 
observe that it is the inclusion of the Gaussian pref actor exp(— a 2 q 2 /2) in (^) which decreases 
the statistical fluctuations dramatically. A small bin width e, which corresponds to a large value 
of a, leads to large fluctuations of (plf), but in the "quantum" algorithm the Gaussian prefactor 
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exp(— a 2 q 2 /2) switches on just in the regime of "small bin width" and thereby dampens out 
the fluctuations. 



4 Two-Particle Bose-Einstein Correlations in a Parton- 
Shower Monte Carlo 

Having gained some insight in the simulation of Bose-Einstein effects within the toy model of 
the previous Section, we now apply the afterburner algorithm to the realistic case of particle 



emission in e + e annihilation at LEP 1 [ |/7j and LEP 2 [ |48fl . We focus on the following reaction 
channels, illustrated in Fig. ||: 

e + e~ -> Z° -> qq -> hadrons at ^fs = 91.5 GeV , (42) 
e + e- -> W + W~ -*qq'q'q -> hadrons at ^=183 GeV. (43) 




Figure 3: Schematics of the two e + e~ event types ( |42| ) and fl4~3"D: The final-state hadron 
distribution in Z° events (left) is due to exclusively 'endogamous' hadronization of the partonic 
offspring from the qq dijet, whereas in W + W~ events (right) there is, in addition, the possibility 
of 'exogamous' hadron production involving a mating of partons from the two different W + — ► 
qq' and W~ — ► g'g dijets. 

These two processes are of interest for several reasons: 

• Generally, e + e~ collisions at \J~(s) > 90 GeV provide the "cleanest" environment of all 
high-energy particle collisions for studying the physics of Bose-Einstein correlations, be- 
cause there is no background to the interesting particles emitted from the calculable par- 
ton shower, and final-state hadrons escape unscathed from their emission point without 
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further interactions. Correlation measurements can therefore be very valuable, as they 
may be used to calibrate analogous analyses in the extreme opposite case of heavy-ion 
collisions, where the emission region is more difficult to calculate accurately, and final- 
state hadron scattering and cascading can crucially influence the shape of the particle 
distributions. 

The experimental study of e + e _ —>■ Z° —>■ hadrons at LEP 1 is based on several million 
events, and hence is impressively extensive and accurate. In particular, high-precision 
measurements of two- and three-particle correlations have been reported [ On the 
other hand, the reaction e + e~ — » W + W~ — > hadrons is currently under very active study 



at LEP 2, in both its experimental and theoretical aspects: for an overview, see [ |48[. The 
interest in this reaction stems from its importance for measuring the triple-gauge-boson 
couplings and mw- in particular, it has been argued that Bose-Einstein correlations may 
introduce an important source of systematic error into the analysis of mw- 

Theoretical studies of Bose-Einstein enhancements have mainly been within the context 
of the string models [ p]], which have been very successful in explaining the distributions of 
identical particles seen in high-energy e + e~ collisions [ ITTJ. Although the string description [ 
[I|] of the hadronization process is a very appealing phenomenological approach and also has 
many other successes, it is not the only possible description. We employ a rather distinct 
cluster hadronization model, based on a space-time description of the perturbative develop- 
ment of parton showers, combined with a non-perturbative model for cluster formation and 
hadronization [ Q . 

The crucial physics point is that, whatever model one uses for the details of the conversion 
of colored partons into color-neutral hadronic states, the Bose-Einstein correlations measured 
in e + e~ experiments are sensitive to local volumes of the order of a fermi in both the longitu- 
dinal and transverse directions. Therefore they provide important information on the intimate 
space-time structure of the hadronization mechanism. In particular, the sources that emit the 
final-state pions and other particles must be identified with local hadronization 'patches', and 
not with the system as a whole, which may extend over even hundreds of fermi. In the string 
picture, these local patches are the centers of string fragments, whereas in our cluster description 
the patches are elementary color-neutral clusters formed from the mating of nearest-neighbor 
partons. The effective Bose-Einstein correlation length should correpond to the sizes of these 
patches, namely the typical string extension ~ 1 fm or the mean cluster size ~ 0.8 fm. Loosely 
speaking, this correlation length defines the minimum possible distance that one may resolve 
from the particle distributions of the hadronic final state. Before turning to our model-specific 
analysis of the Bose-Einstein effect in e + e~ collisions, we refer the interested reader to the com- 



prehensive overview [ |49fl , in which the status of related experimental and theoretical research 
in e + e~ physics can be found. 
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4.1 Modelling the space-time development of e + e collisions 

In order to analyze the effects of identical-particle correlations in e + e _ collisions using the 
quantum version of the Bose-Einstein afterburner, we need to concentrate on event generators 
that deliver realistic, though classical, phase-space distributions of final-state hadrons, to which 
we may then apply the afterburner simulation of quantum interference and the Bose-Einstein 
effect. It is clear from the preceding Sections that such an event generator must not only give 
the momentum spectra, but also the vital space-time information on the dynamical evolution 
and in particular on the final stage of hadron emission. Unfortunately, most of the advanced 



event generators in particle physics [ |5(J do not encode the relevant particle emission structure 
in space and time, whereas most event genarators for heavy-ion collisions do, but cannot be 
applied to e + e _ physics. One event generator that does satisfy both these requirements is 
VNI [[18[], which simulates the e + e~ collision dynamics all the way from the hard annihilation 



vertex, through the perturbative QCD shower development to the emergence of hadronic final 



states. Within the framework of relativistic quantum kinetics | plj, the event generation in VNI 
traces in both space-time and momentum space the parton-shower evolution from the initial 
quark-antiquark pairs, followed by the clustering of the emitted quark and gluon offspring to 
pre-hadronic cluster states that then decay into the final-state hadrons. Referring to [ |, [18|] 
for details, we recall briefly here the essential concepts of this space-time model: 

(i) The parton-shower dynamics is described by conventional perturbative QCD evolution 
Monte Carlo methods, with the added feature that we keep track of the spatial de- 
velopment in a series of small time increments. Our procedure implements perturbative 
QCD transport theory in a manner consistent with the appropriate quantum-mechanical 
uncertainty principle, incorporating parton branching due to real and virtual quantum 
corrections involving gluons or quark-antiquark pairs. In the rest frame of the Z° (for 
(|42|)) or of the (for (fEf)), each off-shell parton i in the shower propagates for a time 
Ati given in the mean by < Ati >= ^iTi = Ei/kf = XiM/2kf, where kf is the parton's 
squared-momentum virtuality, and = E^/M is its longitudinal energy fraction, during 
which it travels a distance Ar^ = Atifii, where M = M z o for fl42[ ) and M = M\y for 
It has been shown [ [| that such a description results in a typical inside-outside 



perturbative cascade [ 52 



(ii) The parton-hadron conversion is handled using a strictly spatial criterion for confinement, 
with a simple Ansatz for the probability P(r) that nearest-neighbor color charges coalesce 
to color-neutral clusters in accord with their color and flavor degrees of freedom, where 
r is the relative distance in between them in their center-of-mass frame. The nearest- 
neighbor criterion is imposed at each time step in the shower development, in such a 
way that every parton that is further from its neighbors than a certain critical distance 
R c = 0.8 fm has a probability distribution smeared around R c for combining with its 
nearest-neighbour parton to form a pre-hadronic cluster, possibly accompanied by one or 
more partons to take correct account of the colour flow. It is important to stress that 
at no moment in this shower development do we incorporate any prejudice regarding the 
genealogical origin of the partons: an 'exogamous' [ |10[ pair of partons from different 
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mother qq' pairs have the same probability of coalescing into a hadronic cluster as do an 
'endogamous' pair of partons from the same mother, at the same spatial separation. The 
resulting hadronic clusters are then allowed to decay into stable hadrons according to the 
particle data tables. 

In the present paper, due to the less demanding requirements on event statistics, we will use 
the "quantum" version of the Bose-Einstein afterburner. This requires fixing the wave packet 
width a in the algorithm. Lacking convincing arguments for a unique physical choice of this 
parameter, we try to connect it with the intrinsic size R c of the pre-hadronic clusters which act 
as pion-emitting sources. The size R c also defines the minimum distance of adjacent clusters 
at formation without overlapping. We thus set 

R c = 0.8 fm = a . (44) 

Finally, we stress that, although the event generation of e + e~ collisions along the above lines 
should provide a rather realistic simulation of the particle dynamics, we do not claim our 
results to be more than qualitative at this point, mainly since we have not included final-state 
interactions among the produced hadrons due to either Coulomb or strong interactions. In 
our approximation, the production vertex of each final-state hadron marks the last point of 
interaction, beyond which the particles stream freely on classical trajectories. Since in 'real- 
life' experiments these final-state interactions can become large at small relative momenta, one 
should be careful when comparing to measured data, some of which have been corrected for 
final-state interactions, others not. 

4.2 Results for two-pion correlations 
4.2.1 Multiplicity distributions 

It is plausible that the structure of the hadronic final state in e + e~ — > W + W~ — > qq' q'q may 
not merely be a copy of W ± — > qq' with twice the final-state multiplicity. As discussed in detail 
in [ [l(|, |{|, it was found within our space-time parton-shower model that not only the total 
multiplicity N(W + W~) may be smaller than 2 x W r± ) , but also that the particle spectra 
may exhibit characteristic differences. These differences are due to the special geometry of 
W + W~ events, in which the partonic offspring of the W + dijet overlap in space-time with the 
partons emitted from the W~ dijet. The cross-talk between the quanta from the is espe- 
cially prominent at small rapidities and if the two dijets emerge at small relative angles. Then, 
whereas in e + e~ — > Z° decays all particles come from the same mother and only 'endogamous' 
cluster formation is possible, as in the left part of Fig. |, W + W~ events receive a significant con- 
tribution from the coalescence of partons from different W mothers into 'exogamous' clusters, 
as in the right part of Fig. 

Fig. [| reflects the effects of parton 'exogamy' in e + e~ — > W + W~ — > qq' q'q as compared to 
Z° — > qq, in the multiplicity distributions of both single pions (top) and of pairs of identical 
pions (bottom). Even allowing for the slightly larger mass of the Z° compared to the W^, 
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Figure 4: Multiplicity distributions of single pions (top) and of pairs of identical pions (bot- 
tom), per charged pion species. 

one observes, in agreement with the above discussion, an effective reduction N 7T ±(W + W~) < 
2xN n ±(Z°), namely, (N W ±(Z )) = 6.1 versus (N 7T ±(W + W^)) = 9.5, per pion species. Similarly, 
we find for the number of identical-pion pairs N w ± n ± = N n ±(N w ± — 1) that N n ± 7T ±(W + W^) < 
4 x N n ± w ±{Z°), namely, (N n ± w ±(Z )) = 81.9 versus (N^±{W + W-)) = 207.3. 

The effect may be thought of as reflecting increased 'efficiency' in the hadronization process, 
due to the fact that the presence of two cross-talking dijets in the W + W~ decays, with their 
spatially-overlapping offspring, allows the evolving particle system to reorganize itself more 
favorably in the cluster-hadronization process, and to form clusters with smaller invariant mass 
than in the Z° events. Indeed, it was found in [ [HJ that the mass spectrum of pre-hadronic 



clusters from coalescing partons is in fact softer in the W + W~ case, reflecting the fact that the 
availability of more partons enables clusters to form from configurations with lower invariant 
mass than in the Z° case. 
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4.2.2 Origins of pions 



In theory, all pairs of identical pions can exhibit Bose-Einstein correlations. Experimentally, 
however, the measurements of the pair spectrum in the relative pair momentum q run out of 
statistics because the phase space vanishes at very low q. Since small q values correspond to 
large spatial distances, this region of q is particularly sensitive to the decays into pions of long- 
living resonances, and also to long-range Coulomb or strong final-state interactions among the 
particles. Whereas final-state interaction effects can be corrected, this is not easy for resonance 
decays. Since many of the pions in e + e _ collisions have their origins in the decays of other 
particles with lifetimes significantly greater than a few fermis, it is useful to disentangle the 
various experimental sources of pions (or other particles) and to classify their parents as follows [ 

• Prompt production leading to pions that emerge directly from the hadronization of the 
fragmenting system, whose parents may be visualized as decaying strings or (in our case) 
as pre-hadronic clusters. 

• Short-lived particles such as p, K* and A, that are strongly-decaying particles with decay 
lengths shorter than a few fermis. 

• Long-lived resonances, such as r\, rj', u, 0, that are states which also decay strongly but 
have life-times of many fermis. 

• 'Stable' particles, such as A and K®, that are particles which propagate sufficiently far 
that the pions emerging can be removed by track cuts. 

• Weakly- decaying particles, such as charm or bottom mesons. 



origin 


life-time r fraction 


clusters 
p,A,K* 


< 0.5 fm 0.31 
1.3 - 4 fm 0.41 
> 10 fm 0.28 



Table 1: Relative contributions of different sources of pions in our e + e event simulations. 

In Table [l], we list the fractions of pions coming from these different sources, as estimated 
in our model simulations. Since the 'stable' particles can be considered as not having decayed, 
and weakly-decaying particles contribute only a negligible fraction, we do not include these two 
categories in the list. We observe that the numbers in Table [I] are very similar to those reported 
in \M 
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4.2.3 The pion correlator C(q, K) for K = 



Fig. [5] shows the correlator C (q, K) — 1 for different q- values and vanishing pair momentum 
K in the cm. frame of the collision, C(q z , q s , q , 0) — 1. Two interesting observations can be 
made immediately: 
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Figure 5: The correlation function of same-sign pions for different values of the relative pair 
momentum q for vanishing pair momentum K, C(q z , q s , q a , 0) — 1. 



a) In both cases, e + e~ — > Z° — ► qq — > n's and e + e _ — > W + W~ — > qq' q'q — > n's, we 
see no significant differences between the three relative momentum directions q z ,qs,<lo- 
Although, for a fixed direction % the intercept of the correlator at qi = depends on the 
magnitude of the momentum transverse to qi, it looks the same for qi = q z , qi = q s or 
qi = q a - 
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b) The correlation function in the case e + e~ — > W + W~ — > gg' g'g — > 7r's is slightly narrower 
than that of e + e~ —>■ Z° —>■ qq —>■ it's. Since the mean q values correspond to the inverses 
of the typical emission source sizes, this means that the hadronic W + W~ decays reflect 
a larger source size (~ 1 fm) than the Z° decays (~ 0.8 fm), as we discuss later in the 
context of resonance effects. 

The implication of observation a) is that the pion emission appears essentially spherically 
symmetric with respect to the three orthogonal directions q z ,q s ,q . This may appear to con- 
flict with the naive expectation that the source should appear much more elongated in the 
longitudinal z direction than in the sideward and outward s, o directions, because of the large 
longitudinal momenta of the leading quark jets. However, as we pointed out before, in this 
model the pre-hadronic cluster formation is controlled by the 'nearest-neighbor' criterion, so 
that only spatially adjacent partons with a mean separation R c ~ 0.8 fm have a significant 
probability of coalescing and decaying into pions and other hadrons. This local coalescence 
results naturally in a longitudinal-momentum ordering of particles as a function of their dis- 
tance from the jet origin: particles further away tend to have higher momentum than those in 
the center. Since the Bose-Einstein effect is only apparent for identical particles with similar 
momenta, corresponding to small q, particles that are separated by many fermi at production 
are incapable of showing a significant enhancement because their momenta are so different. 



One may conclude from observation b) that parton 'exogamy' in W + W~ decays [ 10] results 
in a space-time distribution of hadrons that is more spread out than in the case of Z° decays. 
This may again be understood as a consequence of increased efficiency of hadron formation in 
the W + W~ events, as we discussed before in the context of the pion multiplicity distributions in 
Fig. |J The identical pions emerging as products of parton-cluster decays have a longer distance 
correlation in W + W~ events, because the partonic offspring of the overlapping qq dijets are 
enhanced mainly for low-momentum quanta in the central rapidity region, corresponding to 
significant 'exogamous' coalescence of partons from different W's with small relative momenta. 
As a consequence, the pion pair spectrum from W + W~ decays in Fig. [|b is narrower than the 
one from Z° decays in Fig. ||a, which translates into a larger effective emission radius for these 
pions. 



4.2.4 The pion correlator C(q,K) for K ^ 



The general features and physics interpretation of the Independence of the correlation function 
C(q,K) have been studied in detail in [ 25, |SJ| |53{] . A manifest change in the shape of 



C(q,K) as K varies can have several origins, the two most important being (i) resonance 
decay contributions [ ^ and (ii) collective flow of the particle matter [ [ZBj. Whereas pions 
from long-lived resonances are always present in high-energy collisions, collective motion of 
the produced particles is a feature of heavy-ion reactions that produce high-density matter, 
but certainly is not an issue for the e + e~ collisions discussed here. In this subsection we do 
not distinguish the contributions from resonance decays, but show the Independence of the 
correlation function including pions from long-lived resonances, just as in the previous figures. 
We disentangle the effect of resonance decays in the next subsection. 
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Fig. |^ shows the correlation function C(q,K) — 1 of same-sign pions for various values of 
the pair momentum K = (Kl,K±), where Kl = K z is the direction along the thrust axis 
and K± = \K±\ = JK% + K$ the momentum transverse to it. The correlator is plotted as a 
function of one of the three Cartesian components of the relative momentum q, with the other 
two components set to zero. The two main features are: 
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Figure 6: The correlation function of same-sign pions C(q,K) — 1 for various values of the 
pair momentum K = (Kl,K±), where Kl = K z is the direction along the thrust axis and 
K± = \K±\ = \JK^ + Ky is the momentum transverse to it. The correlators are plotted against 
one component of the relative momentum, setting the two other components to zero. 



a) The shape of the correlation function flattens and widens as Kl or K± are increased. 
However, the mean q^x change by less than 10 % when the Kl,t values are varied from 
to 1 GeV. 
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b) The /^-dependence is evidently spherically symmetric, i.e., the correlation function changes 
in the same way as Kl or K± is increased in the range from to 1 GeV. 



Point a) is a reflection of resonance-decay pions: long-lived resonances with r > 10 fm can 
travel for many fermis before decaying, which leads to an exponential tail in the pion pair 
spectrum [ This 'life-time effect' is larger for small values of K and damps out as K is 



increased, since the relative abundance of resonances is most pronounced at small K. 

Point b), on the other hand, is in accord with the fact that the kinematics is approximately 
boost-invariant along the thrust axis of the e + e~ collision, and the 'local' character of the 
particle dynamics in our model. Neither the parton shower evolution nor the parton-cluster 
hadronization depend on the overall momentum K relative to the e + e~ center-of-mass frame, 
as it is only the kinematics, color and flavor of the near-by clustering partons which at any 
given vertex determines locally the development of particle production. 



4.2.5 Effects of resonance decays on C(q,K) 

Consider now a pair of identical pions with relative momentum q, where one of the pions 
originates from a resonance of momentum p with mass m r and decay width T r ~ r" 1 . Such 
a pair cannot contribute to the Bose-Einstein effect if \q ■ p\ 3> m r T r , which roughly implies 
that \q\ 3> r" 1 . Since q is inversely proportional to the spatial dimension of the pion source, 
this means that resonances represent a source of spatial extent of the order of r r _1 . Hence, such 
pions only can exhibit correlations if \q\ < 0(r r _1 ), which for long-living resonances (r r > 10 
fm) requires \q\ < 20 MeV. This is less than the scale at which direct pions coming from the 
pre-hadronic clusters contribute, as seen in Table [I]. Therefore the pion correlation function 
C(q, K) is narrowed by the effects of the pion decay products of long-living resonances. 

To quantify the resonance narrowing and localize the pile-up of pion pairs where at least 
one comes from resonance decay, we have disentangled the pion emission sources of Table p] 
within our event simulation. In Fig. [7] we show again the correlator C(q, K) — 1 for K = in 
the two cases of e + e~ — > W + W~ — > qq' q'q — > 7r's and e + e _ — > Z° — > qq — > 7r's. The full lines 
correspond to the correlator with all sources included, as in the previous Fig. [], whereas the 
dashed curves have the long-living resonance decay contributions removed. One sees that the 
resonance decay pions make a significant 20 - 30 % contribution to the magnitude of C — 1 at 
Ql-iQt ~ 50 MeV, corresponding to life-times > 5 fm. Fig. |8| exhibits the Independence of the 
effective emission radii associated with the mean values of the components of q = (g^, q s , q ) in 
the longitudinal (Rl) and transverse (R ,R S ) directions with respect to the thrust axis [ ^4j: 

Again, the full lines include all the sources in Table |l|, whereas the dashed lines exclude the 
long-living resonances. It is evident that the effect of resonances is most pronounced for small 
values of K and disappears with increasing K. This is expected as the abundance of pions 
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Figure 7: The correlator C(q,K) — 1 for K = with (solid curves) and without (dashed 
lines) the contributions from long-lived resonances. 

from resonance decays is most prominent at small K. Resonance decays thus induce a pair- 



momentum dependence of the HBT radii [ 0, [55], with overall variations of the radii on the 
order of 0.1 fm. 

It is interesting to note that after switching off the resonance decay contributions the HBT 
radii shown in Fig. [8] exhibit no remaining pair- momentum dependence. Such a pair- momentum 
dependence would signal the presence of position- momentum correlations in the source as, e.g., 
induced by collective flow or string breaking kinematics. No such correlations are visible here, 
not even along the thrust axis (see the left panels of Fig. [| which show Rl as a function of K±). 
This is surprising because the inside-outside cascade features of parton and hadron production 
in VNI should lead to appreciable position-momentum correlations along the longitudinal axis 
defined by the primary hard partons. We can think of two possible reasons for the fact that they 
are not reflected in the longitudinal HBT radii shown here: either they get largely averaged 
out by summing over many collision events (which we think is unlikely), or they get covered 
up by the finite size of the wave packet width a in the Bose-Einstein afterburner. The fact 
that all HBT radii come out very close to o = 0.8 fm lends support to the second conjecture, 
although a final clarification of this issue has to await a comparison with calculations based 
on the "classical" version of the afterburner, as well as studies of the "quantum" version with 
different values of a . 
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Figure 8: The K dependence of the effective pion source radii associated with the mean values 
of the components of q = (g^, q s , q Q ) in the longitudinal {Rl) and transverse (R Q , R s ) directions 
with respect to the thrust axis. The solid lines include all sources, whereas the dashed curves 
exclude long-lived resonances. 



5 Discussion 



We have discussed in this paper two possible algorithms for modelling Bose-Einstein correlations 
in a Monte Carlo code for e + e~ annihilation into hadrons, that incorporates information from 
perturbative QCD on the space-time evolution of parton showers and a configuration-space 
criterion for hadronization. Afterburners incorporating both the "classical" and "quantum" 
algorithms have been applied to a model in which the hadron emission region is known ana- 
lytically. Standard tools for analyzing the sizes of hadron emission regions have been applied 
to these model calculations, and shown to reproduce successfully the parameters of the model. 
The quantum algorithm has then been implemented as an afterburner in the space-time parton- 
shower Monte Carlo, and applied to e + e~ — > Z° — > hadrons and e + e~ — > W + W~ — > hadrons. 
Exploratory analyses have been presented of two-pion correlations in longitudinal and trans- 
verse momenta, both with and without resonance decays. The latter have been shown to modify 
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significantly the Gaussian behaviour that would otherwise have been expected, and to cause a 
pair- momentum dependence of the extracted HBT radii, albeit on a small scale of order 0.1 fm 
only. In the limit studied here where the "quantum" afterburner was used with a fixed wave 
packet width a = 0.8 fm, resonance decays in fact induced the only discernible K±- dependence 
of the HBT radii. 

The analysis of this paper has necessarily been incomplete, and we conclude by listing some 
of the open questions that could be addressed in any future work. It would be interesting 
to implement the "classical" algorithm as an afterburner, and investigate the similarities and 
differences with the quantum afterburner explored in this paper. We are not in a position 
to express a definitive theoretical preference for one algorithm over the other. Within the 
context of the quantum afterburner, we have assumed one particular value of the Gaussian 
wave-packet size a, and have not explored the implications of varying this parameter. In this 
connection, we should draw the reader's attention to the possibility that the similarities between 
the correlations in the transverse and longitudinal momenta may be related to the choice of a, 
which we have not attempted to optimize. This would involve an overall tuning of the Monte 
Carlo to fit particle spectra, which we are currently not in a position to complete^. Once this 
is done, one could use the Monte Carlo to address some of the physics issues that triggered this 
investigation, including the possible effects of Bose-Einstein correlations on measurements of 
the mass in hadronic final states at LEP 2. 

Despite the inevitable incompleteness of this work, we hope that the ideas and investigations 
reported here may be useful in future studies along the lines suggested above, either within 
the context of the space-time approach to e + e~ annihilation into hadrons used here, or within 
some other approach. The Bose-Einstein afterburners studied here could also be implemented 
in Monte Carlo codes for other interactions, including relativistic heavy-ion collisions. We are 
currently studying how this work could be advanced along these lines. 

The authors acknowledge illuminating discussions with M. Gyulassy, J. Sollfrank, S. Vance 
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